Thermal conductivity measurements of porous dust aggregates: 
I. Technique, model and first results 



M. Krause^ J. Blum^ Yu. V. Skorov^ M. Trieloff* 

^Institutfur Geophysik und extraterrestrische Physik, Technische Universitdt Braunschweig, 
Mendelssohnstr. 3, D-3 8 106 Braunschweig, Germany 
^ Institut flir Geowissenschaften, Ruprecht-Karls-Universitdt Heidelberg, 
Im Neuenheimer Feld 234-236, D-69120 Heidelberg, Germany 



o ■ 

Abstract 

M We present a non-invasive technique for measuring the thermal conductivity of fragile and sensitive mate- 
^ rials. In the context of planet formation research, the investigation of the thermal conductivity of porous 
^ dust aggregates provide important knowledge about the influence of heating processes, like internal heating 
I— ^by radioactive decay of short-lived nuclei, e.g. ^^Al, on the evolution and growth of planetesimals. The 
determination of the thermal conductivity was performed by a combination of laboratory experiments and 
numerical simulations. An IR camera measured the temperature distribution of the sample surface heated 
^ by a well-characterized laser beam. The thermal conductivity as free parameter in the model calculations, 
^ exactly emulating the experiment, was varied until the experimental and numerical temperature distributions 
^ showed best agreement. Thus, we determined for three types of porous dust samples, consisting of spherical, 
^ .1.5 yum-sized Si02 particles, with volume filling factors in the range of 15% to 54%, the thermal conductivity 
I — ^!to be 0.002 to 0.02 Wm"^ K"\ respectively. From our results, we can conclude that the thermal conductivity 
^ mainly depends on the volume filling factor. Further investigations, which are planned for diff'erent materials 
> and varied contact area sizes (produced by sintering), will prove the appropriate dependencies in more detail. 

OO Keywords: Comets, dust. Comets, origin. Origin, Solar System, Planetary Formation, Planetesimals, 
Regohths 
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^ 1. Introduction 
^ ! 

I In the field of planet-formation research, the in- 
vestigation of the physical properties of porous dust 
^ aggregates - as analogs for protoplanetary bodies 
. 5^ - is of great importance for understanding the evo- 
lution of solid matter in young planetary systems. 
The evolution of the size (distribution) of dust aggre- 
gates depends heav ily on the collision properties of 
these p articles (see Zsom et al. ( 2010b : Guttler et al. 
(|2010l) \ For larger dust aggregates, other proper- 
ties become also relevant. Here, we will discuss 
the thermal conductivity of macroscopic dust aggre- 
gates with various packing densities and try to fig- 



Email address: maya . krause@tu-bs . de (M. Krause) 
Preprint submitted to Icarus 



ure out to what extent heating processes can alter 
their internal structure. Beside globally elevated tem- 
peratures in protoplanetary disks and episodic heat- 
ing events, internal heating of larger bodies by ra- 
dioactive decay of short-lived nuclei, e.g. ^^Al, can 
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Planetesimals are initially highly porous and frag- 
ile bodies. As high porosity bodies have only few 
contacts between neighboring monomers of the ma- 
terial and/or very small contact areas between them, 
a lower thermal conductivity is expected than in 
more compact bodies. If a porous body is heated 
from inside by radioactive decay, sintering between 
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monomer grains in contact can occur. By the for- 
mation of inter-particle necks due to the sintering 
process, the contact areas between the single parti- 
cles will be enlar ged and thus the material solidi- 
fies (IPoppd . I2OO3I) . The larger contact areas of the 
sintered material would imply a higher thermal con- 
ductivity and as a consequence a more efficient and 
faster heat transport to the surface of the body. The 
intimate interplay between heating due to low ther- 
mal conductivity and an increase in heat flow caused 
by sintering is a yet unsolved problem in planet- 
formation research. The question whether these two 
counteracting eff'ects lead to a complete or partial 
melting of the body during growth and how such ef- 
fects influence the further growth and the structural 
and physico-chemical evolution of planetesimals are 
yet to be answered. 

Although the thermal conductivity of porous 
media is an important quantity for heat trans- 
fer in many other di sciplines, like e.g. g eolog y 
and engineering (e.g. van Antwerpen et al] ( 2010 ): 



Masamune and SmithI (| 19631) ). comparable thermal 



conductivity measurements for highly porous mat- 
ter do not exist, which vary the bulk porosity and 
packing structure, use particle sizes in the order of 
few micrometers, and are performed under vacuum 
conditions. Most of the literature about heat con- 
duction of particulate materials deals with gaseous 
or fluid transport through the voids, whereas typical 
gas pressures in protoplanetary disks or (at a later 
stage) debris disks are so low that the heat trans- 
port through the pore space of dust aggregates by 
gas molecules can be neglected (see also Sect. (6]). 
Comparable measurements under vacuum conditions 
or at low atmospheric pressures found in the litera- 
ture mostly concentrate on larger particle sizes and 
higher (a nd often inhomogen e ous) packing densi - 



ties (e.g. iHuetter et all (120081) : iKtihrt et all (119951) : 



Merrilll (119691) ). 

The presumably most unaltered protoplanetary 
material available in the Solar System can be 
found in comet nuclei, and recent infrared mea- 
surements of Comet 9P/Tempel 1 have shown that 
the thermal inertia / = ^jkpC, with k, p, and 
C being the heat conductivity, the density, and 
the specific heat capacity, r espectively, is as sma ll 
as / < 50WK-^m-2s^/2 dOroussin et al.L l2007b . 



Assuming p = 100 . . . 1, OOOkgm"^ and C = 
700...1,400Jkg-^K-\ a thermal conductivity of 
< 1.8 X 10-^..3.6 X lO-^Wm-^K-i can be 
derived. Mind, that the low thermal-inertia val- 
ues derived by lOroussin et al.l (120071) are in con- 
flict with a diff e rent a nalysis of the same data by 



Davidsson et al.l (|2009h . who yield thermal inertias 



of wide surface ranges of Comet 9P/Tempel 1 of 
1000-3000 W m"^ s^/^ whereas only small frac- 



tions of the comet surface possess low thermal iner- 
tias of 40 - 380 W K"^ m"^ s^^^. The measurement of 
the thermal conductivity of porous dusty media can, 
thus, provide information about the morphology and 
structure of the surface material of comets. Such re- 
sults are not only important for the early evolution of 
small Solar System bodies, but also for thermal pro- 
cesses in evolved systems, e.g. for heating and im- 
plicit activity of comets with their porous ice struc- 
ture and dust mantles during passage through the in- 
ner Solar System, or the Yarkovsky eff'ect for aster- 
oids that may have regolith on their surfaces. 

The thermal conductivity of porous media is af- 
fected by a large range of parameters: the bulk 
material in general (pure or mixed), the tempera- 
ture, the porosity, the shape of monomers, the size 
distribution of the grains (mono- or poly disperse), 
the number of contact points between the parti- 
cles, the size distribution of the cont act area (poten- 



tially c h anged by sinter i ng eff ects, see lKossacki et al. 



(119941) : ISeiferhn et all (119951) ). and the structure of 
the bulk matrix. The latter can be disordered or or- 
dered, and in the ordered case can additionally be 
divided into isotropic (with simple or complex struc- 
ture) and anisotropic (isotropic in one plane or not 
isotropic in any plane). 

Many models to describe the evolution of 
planetesimals or parent bodies of meteorites use 
thermal conductivity val ues derived by mea- 
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terial has already undergone several thermally and 
mechanically induced structural modifications and, 
thus, cannot correctly represent the primary highly 
porous state of these bodies. Even an often used 
correction or reduction factor for the heat conduc- 



tivity to account for the porosity or the cohesive- 
ness of the bulk material cannot describe the com- 
plex correlation between the thermal conductivity 
and the i nner structure of the material to the full ex - 
tent fe.g.lCapriaetalJ (l2002b : ISeiferlin et all (Il996b : 



Paton et alJ (|201Qh y For many decades, extensive 
studies have been made to describe the relationship 
between the thermal conductivity and the porosity 
of two-p hase mixtures and porous materials (see re- 
views by IProgelhof et alJ (119761) . ICheng and Vachon 
(I197Q|) \ All the semi-empirical equations mentioned 
in the literature are based on well-defined mate- 
rial configurations and environmental conditions and 
thus are restricted to specific application problems. 
In addition to that, most of these formulae consider 
the porosity as the most important parameter defining 
the thermal conductivity without including the influ- 
ence of pore size, pore shape, material structure, etc.. 
In this work, the numerical model reproducing our 
experimental measurements does not include a de- 
tailed theoretical description about the thermal con- 
duction of porous material but uses a reduction factor 
of the thermal conductivity as a free model parameter 
(see Sect. |4l), which indirectly but not explicitly in- 
cludes the manifold structural parameters of the bulk 
material contributing to the heat conduction. 

To understand the intricate dependency between 
the thermal conductivity and the structure of the pro- 
toplanetary body, extensive experimental investiga- 
tions are needed. The results of these measurements 
will serve as important material parameters for mod- 
eling the internal constitution and thermal evolution 
of planetesimals and cometary nuclei and generally 
for other thermal proce sses such as thermophor esis 
and photophoresis (e.g. IWurm and Haac^ toO^ ). 

For investigating the heat conductivity of porous 
and fragile dust samples, it is very important to 
choose an appropriate measuring technique. Con- 
ventional methods, like the hot wire or the hot 



plate m ethod (see review by lPresley and Christensen 



all need direct contact to the material. By 
establishing the contact between the measuring appa- 
ratus and the material of loose dust assemblies, the 
frail inner structure of the dust sample will always 
be changed in terms of compression or disruption. 
Compressive and tensile strengths of our materials 
are typically in the range of a few 100 to a few 1,000 



Pa (IBlum et al.L l2006r) . and even the slightest con- 
tact with a heating wire will cause unknown local 
changes of the aggregate morphology. Thus, we have 
developed a non-invasive technique to measure the 
thermal conductivity by heating the dust aggregate 
with a laser beam and recording the temporal and 
spatial propagation of the heat wave by an IR camera, 
which is comparab l e to the laser flash tec hnique (e.g. 
Parker et all (Il96lb : iFayette et all tood) ). The tem- 
perature distribution on the surface and in the interior 
of the sample is modeled using variable thermal con- 
ductivity values until the best fit to the experimental 
results is achieved. 

In this first publication about thermal conductivity 
measurements in the context of planet formation re- 
search, we will present the measuring technique and 
the corresponding model calculations for a set of ex- 
emplary dust samples with diff'erent porosities and 
diff'erent internal structures, consisting of the same 
dust type. Further measurements with other pro- 
toplanetary relevant materials and diverse sintering 
stages are under way. 

2. Experimental Setup 

As we expect very low values for the thermal con- 
ductivity of porous dust aggregates in the order of 
~ 10"^ . . . 10"^ Wm"^ K"\ as indicat ed by the mea- 
surem ents of Comet 9P/Tempel 1 (IGroussin et al 



(|2007|) . see above), high- vacuum conditions are re- 
quired for the measurements, as the presence of 
(even rarefied) air would significantly contribute to 
the thermal conductivity. Hence, all measurements 
were performed in a vacuum chamber at pressures 
around 10"^ mbar. From outside the vacuum cham- 
ber, an infrared laser beam (wavelength = 813 nm), 
an overview camera to align the laser, and an IR cam- 
era are pointed through respective windows onto the 
aggregate inside the chamber. To prevent additional 
sources for heat conduction, the dust sample is posi- 
tioned onto a block of insulation material with a ther- 
mal conductivity of ~ 0.022 Wm"^ K'^ The thick- 
ness of the dust samples is so large that we do not 
expect a considerable heat flow to this substrate (see 
Sect. [5] and Fig. [3]). The cylindrically shaped dust 
aggregate (see Sect. [3]) is directly heated by the laser 
beam mounted vertically to the surface of the dust 



sample (see Fig.[T]). 
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Figure 1: Experimental principle for measuring the thermal 
conductivity of a cylindrically shaped dust sample. The sam- 
ple is heated by a laser beam, while an IR camera monitors the 
temporal and spatial temperature distribution. 



well-known to us in terms of physical properties after 



years of research (|Blum et all 120061 : iGtittler et al. 



2OIOI) . One additional advantage to choose this mate 



rial for our initial experiments is the spherical shape 
of its monomers and their well-defined size. For fu- 
ture numerical simulations, information concerning 
the internal structure of the sample material is re- 
quired to model the heat conduction through the con- 
tact areas of neighboring individual particles. The 
geometrically simple shape and structure of this dust 
material facilitates the development and testing of 
such a numerical thermal conductivity model. 

To study the influence of the material porosity 
on the thermal conductivity, we manufactured and 
analyzed dust samples of diff'erent porosities, re- 
sulting from three diff'erent production techniques. 
By the random ballistic deposition (RED) method 



During the heating and cooling phase, an IR cam- 
era monitors the temporal and spatial temperature 
distribution of the aggregate's surface at an inclina- 
tion angle of ~ 30° to the surface normal. The im- 
ages of the IR camera have a spatial resolution of 
~ 0.3 mm and a temporal resolution of 0.02 s. The 
temperature measurements of the IR camera were 
performed with a default emissivity value of unity. 
To achieve more accurate temperature values, we cal- 
ibrated the temperatures of the IR camera by compar- 
ison to a heated material with the same IR emissivity 
value as the dust sample, using a thermocouple for 
precise direct temperature measurement. Each sam- 
ple was heated by the laser beam until the tempera- 
ture of the spot of the laser beam on the sample sur- 
face did not considerably change, i.e. the heating by 
the laser and the conduction of the heat inside the 
dust sample were reaching an equilibrium state. The 
cooling phase, after the laser had been switched off', 
was recorded until the dust material of the sample 
was cooled down to room temperature. 

3. Sample Preparation 

For the first thermal conductivity measurements, 
we used Si02 dust, consisting of monodisperse 
spheres with a diameter of 1.5 yum. This material 
is widely used as protoplanetary dust analog and is 




Figure 2: Examples of dust aggregates manufactured by the 
three different methods described in the text. Top: dust sam- 
ple produced by the RBD method; middle: sieved dust sample; 
bottom: compressed dust sample. 



(|Blum and Schrapleil . 120041) we produced highly- 
porous cyhndrically shaped dust samples with a di- 
ameter of 2.5 cm and a volume filling factor of = 
0.15 (see Fig. [2l top). To achieve moderate vol- 
ume filling factors in the range of ^ = 0.15 .. . 0.30, 
we sieved the dust through metal filters with diff'er- 
ent mesh apertures into a cylindrically shaped plastic 
container with a diameter of 3.5 cm (see Fig. [2l mid- 
dle). The collection of the sieved dust into the con- 
tainer was necessary to achieve well-defined borders 
of the sample and a flat surface. Previous measure- 
ments showed that the propagation of the heat wave 
for this material is not large enough to reach the bor- 
ders of the container so that the container material 
has no influence on the thermal conductivity of the 
dust sample. More compact samples with a volume 
filling factor of ^ 0.5 were produced by compress- 
ing the dust by a plunger inside a hollow cylinder 
with an inner diameter of 2.5 cm (see Fig.[2l bottom). 



4. Modeling 



The problem to derive the thermal conductivity 
of a body from measured distributions of surface 
temperatures can be solved by numerical modeling. 
The heated sample is described as a straight cylin- 
der with height and radius measured in the experi- 
ment, i.e. the solution to the heat conduction equa- 
tion is found in a two-dimensional axisymmetric ge- 
ometry. The volume filling factor is defined exper- 
imentally, whereas the heat capacity and the den- 
sity of the solid mate rial are taken from the literature 



to be 840Jkg-^K-^ jripL 



ei 



19991) and 2,000 kg m"^ 
(|Blum and Schraplerl . 120041) . respectively. It is as- 



sumed that the sample does not contain volatile com- 
ponents and that the ambient pressure is extremely 
low (because the actual experiment was conducted 
in a vacuum chamber, see Sect. [2]). 

A two-dimensional transient nonlinear heat con- 
duction equation for the porous sample with volume 
filling factor (j) is solved. Note that we use a continu- 
ous description for the porous medium. This means 
that the classical heat transfer equation is used with 
the eff'ective parameters (density p(r) = po0 and 
heat capacity c(r)), defined for a particular porous 



medium 



dt 



c(r)po(r)(f>(ur) = V(/:(r, r)Vr(^, r))+^(^, r) . 

(1) 

The first term on the right-hand side of Eq. [T] de- 
scribes the net change of heat in the unit volume due 
to heat conduction. The eff'ective conductivity k(T, r) 
is given by k(T, r) = h kd+kjR, where the approximate 
bulk conductivity of a solid mixture of dust, ^j, has 
been corrected for the aggregate structure via a re- 
duction factor h. This factor is strongly dependent 
on the morphological structure of the medium, e.g., 
the monomer size and shape, the porosity, the coor- 
dination number, and the sintering state of the dust 
material. As these parameters are hard to express 
quantitatively, we treat as a free parameter in the 
present work. It is assumed that the conductivity of 
solid dust kd is constant and equal to 1.4 Wm"^ K"^ 
We also assume that the dust particles are completely 
opaque in the thermal radiation range. The radiative 
conductivity kjR is taken into account when the to- 
tal thermal conductivity of the medium is calculated, 
and is defined as 



(2) 



with the emissivity s equals (1 -albedo) according to 
Kirchhoff^'s law, cr the Stefan-Boltzmann constant, 
and / the inter-particle distance. 

Since all the samples used in our experiments are 
characterized by low packing densities (the maxi- 
mum value of the filling factor is equal to 0.5 for the 
compressed sample), we include bulk absorption of 
energy in the numerical model. This means that the 
incident radiation is absorbed in the volume, rather 
than on the sample surface. We assume that the 
weakening of the intensity of the direct laser radia- 
tion obeys Beer's law, i.e. decreases exponentially 
with depth. The second term on the right-hand side 
of Eq.[T]is a source term for the laser energy absorbed 
by the unit volume 



^(^r) = £oexp(-r//l) , 



(3) 



here E{t, r) is the net downward directed energy flux, 
E{) and A are a normalization factor and an attenua- 
tion length, respectively. 

We do not consider at this stage the eff'ects associ- 
ated with multiple scattering of direct radiation in the 



porous layer. The fact is that, although our samples 
possess significant porosity, this porosity is not suffi- 
cient to apply the traditional computational methods 
of modeling the transport of electromagnetic radia- 
tion, designed for sparse environments. Application 
of the classical equation of radiative transfer can re- 
sult in noticeable errors already at volume fillin g fac- 
tors of = 0.30 (ITishkovets and JockerslJ2006b . 

The upper boundary condition for the temperature 
is obtained from the energy balance equation for the 
surface. Energy is given to the surface by laser ra- 
diation and heat conduction (in case of positive tem- 
perature gradient). Energy is removed from the sur- 
face due to thermal radiation and heat conductivity 
(in case of negative temperature gradient). Mathe- 
matically, the condition is formulated as 



kVT{t, r) = sctT^ 



(4) 



where we have a balance between thermal re- 
radiation and heat conductivity. We note that the ra- 
diation term is absent in the case of porous media. In 
this way, the heat tra nsfer equation deals with Radia- 
tion as a source term (IDavidsson and SkorovlJ2002r) . 
The thermally insulating boundary condition was ap- 
plied to all the other boundaries of the model. This 
simple condition has been chosen after tested simu- 
lations clearly demonstrated that the heat wave does 
not reach the model boundaries during the experi- 
ment. 

An important part of our numerical simulation is 
the model of the laser beam, which serves as the sole 
source of heating energy. In the computer model, 
the laser beam is characterized by its total power, the 
spatial profile of the intensity, and the duration of 
the illumination. All these characteristics were de- 
termined experimentally (see Sect. 14.11) . 

We now describe the numerical model. Our com- 
puter experiments (as well as our laboratory exper- 
iments) can be divided into two phases: an active 
phase, when the sample is irradiated by the laser 
beam, and a passive stage, after the laser is switched 
off' and the heat (wave) is spreading in the porous 
material. As an initial condition for the temperature 
inside the samples, an isothermal temperature distri- 
bution is used. This initial temperature is given in ac- 
cordance with the specific experiment. The boundary 
condition at the top of the cylindrical sample is the 



energy balance and, under the above assumptions, 
this condition contains only a term that defines the 
thermal radiation of the sample and a thermal flux. 
At all other boundaries of the sample, the insulation 
condition is applied. 

A non-uniform adaptive grid is used in the numer- 
ical simulations. The minimum spatial step is 10 yum. 
Usually the numerical grid contains about 30,000 el- 
ements. We use UMFPACK, which c an solv e linea r 
sparse systems via LU factorization (|Davisl . 120041) . 
This direct solver is applied to the whole matrix of 
grid points and is efficient for our 2D problem, al- 
though it may be too memory-intensive for 3D sim- 
ulations. 

From the point of view of computer modeling, the 
specific feature of the considered problem is the pres- 
ence of three free model parameters, i.e. the heat 
conductivity of the sample, the albedo of the sample 
at the wavelength of the laser beam, and the penetra- 
tion depth of the laser beam into the sample. This 
fact makes the task of determining the eff'ective ther- 
mal conductivity more complex than the standard so- 
lution of the non-linear transient heat transfer prob- 
lem in porous media. In fact, our problem can be for- 
mulated as the problem of finding the optimal set of 
model parameters that yields the best agreement be- 
tween the numerical simulations and the laboratory 
experiment. 

To solve this complex problem, we essentially 
relied on the setup of the laboratory experiment, 
namely its division into an active (heating) and pas- 
sive (cooling) stage. In the first phase, an inhomoge- 
neously heated region is formed in the sample. The 
particular temperature distribution in this zone is de- 
fined by the absorption of the laser energy as well 
as by the heat transfer due to contact and radiative 
conduction. During the cooling stage, the radiation 
source is turned off', and the evolution of the heated 
zone is determined only by the eff'ective thermal con- 
ductivity of the medium. This fact makes the task 
of the derivation of the thermal conductivity easier. 
Thus, initially the problem of the quantitative eval- 
uation of the model parameters Eq^ and A is solved. 
The experimental estimates for the radiation absorp- 
tion and for the size of the area in which almost all 
the radiation is absorbed (see Sect. 14.11) , are used 
as initial approximations in these simulations. The 



temperature evolution during the active phase is char- 
acterized by the first and second time derivatives of 
the temperature and by the maximum surface tem- 
perature achieved at the end of the irradiation phase. 
These characteristics, derived from the laboratory 
and the computer experiments, are compared for a 
few selected values of the distance from the axis of 
symmetry of the sample. Then, using the calculated 
temperature distribution obtained at the end of the 
heating phase as an initial temperature distribution 
for the further simulations, we treat the thermal con- 
ductivity as the only free parameter and make a fine 
adjustment, looking for the best compliance between 
the computer model and the laboratory experiments. 



4.1. Measurements for Model Parameters 

To reproduce the laboratory experiments by the 
numerical simulations in the most realistic way, sev- 
eral additional measurements were necessary. The 
volume filling factors of the diff'erent dust samples 
were determined using the relation (p = m/(Vpo) be- 
tween the mass m and the volu nie of the dust sample 
V. Here, po = 2, 000 kg m"^ (|Blum and Schraplerl . 



20041) is the density of the solid material. In the case 
of the sieved dust samples, the volume filling factor 
of the bulk sample (l)tuik = ^agg (paa comprises the vol- 
ume filling factor of the individual sifted (sub-mm- 
sized) aggregates, (pagg, and the volume filling factor 
of the "aggregate of aggregates", 0^^, which forms 
when the sifted aggregates fall into the repository. 
Since it is very difficult to measure the porosity of 
single very small particles with a diameter in the or- 
der of the mesh apertures a^esh = 0.15 .. . 0.50 mm, 
we estimated these volume filling factors as shown 
in Table [B 



Table 1: Volume filling factors and packing densities of the 
sieved dust samples. 



a„esh [mm] 


<Pbulk 


^agg 


<Paa 




(measured) 


(estimated) 


(estimated) 


0.50 


0.29 


0.45 


0.64 


0.25 


0.24 


0.45 


0.53 


0.15 


0.16 


0.45 


0.35 



As the mesh widths are not too diff'erent in size, we 
assume an identical volume filling factor of (pagg = 



0.45 for the sieved agglomerates. This value re- 
sults if we suppose for the largest sieved agglomer- 
ates a packing density of (paa = 0.64 , which cor- 
responds to a random close packing of equal-sized 
spheres. In fact, the sieved dust aggregates are not 
perfectly spherically shaped and do not exhibit a 
monodisperse size distribution, what would yield a 
somewhat higher packing d ensity than ^aa = 0.64 
(|Schaertl and SillesculJl994|) . In contrast to this, the 
sieving process of dust aggregates leads to a some- 
what less dense packing than the random close pack- 
ing of solid spheres, due to non-negligible cohe- 
sion between neighboring aggregates. As we think 
that these contrary aspects will probably compen- 
sate each other, we decided for a packing density 
of (f>aa = 0.64 for the largest sieved agglomerates 
as an upper limit. In any case, a filling factor of 
(pagg = 0.45 is in agreement with measurements of 
dust aggregates sifted with larger mesh sizes (R. Wei- 
dling, pers. comm.). 

The heating process by the laser beam is charac- 
terized within the simulation model by the total laser 
power, the spatial energy distribution within the laser 
beam, the albedo of the sample, the penetration depth 
of the laser beam, and the duration of illumination, 
which we all empirically derived or estimated. Re- 
garding the laser beam, we measured its profile and 
eff'ective power for diff'erent input voltages. Due to 
an adapted coUimation optics, the laser profile de- 
viates from a gaussian profile at the center and by 
wider flanks; the full width at half maximum power 
of the laser profiles was determined to be ~ 1 .2 mm. 
To determine how deep the laser beam penetrates 
into each dust sample type during the heating, we in- 
vestigated the transmission of the laser light through 
dust layers of diff'erent thicknesses. We found that 
for the RBD samples, a sample thickness of 1 mm 
attenuates the laser intensity to ~ 10"^ of its origi- 
nal value; for the sifted samples, the attenuation after 
penetrating a sample thickness of 1 mm was ~ 10""^ 
and for the compressed samples ~ 10"^. These re- 
sults verify that for all sample types after a depth of 
a few millimeters, absolutely no light is transmitted 
anymore. 

To estimate the amount of laser light absorbed by 
the samples, we made use of specific albedo mea- 
surements performed at the PTB in Braunschweig. 



These albedo measurements were restricted to the 
compressed samples due to the vertical orientation 
of the sample holder and yielded a value of 0.983 
for the albedo. For such high albedo values, how- 
ever, the measurement uncertainties are quite high. 
For the more porous samples, the albedo estimation 
is technically even more challenging. In addition, the 
increase of porosity can lead to significant variations 
in albedo associated primarily with the local inho- 
mogeneity of the surface. Thus, we use the mea- 
sured albedo values for an indication of very high 
albedos and leave the albedo as a free parameter in 
the model simulations. It is again important to note 
that our computer modeling is naturally divided into 
two stages: the stage in which the sample is heated 
by the laser beam (active phase), and the stage when 
the sample cools down (passive stage). In the second 
stage, the free model parameters albedo and absorp- 
tion profile of the radiation inside the sample play no 
direct role. However, the corresponding time evolu- 
tion of the heated sample depends on the volumetric 
distribution of absorbed laser energy and, thus, the 
temperature at the end of the active phase. There- 
fore, albedo and the attenuation profile of the laser 
radiation have an indirect influence on the subse- 
quent thermal evolution of the sample. As will been 
seen in Sect. [51 the values of the free model pa- 
rameters used in our simulations are in the expected 
range of the available experimental estimates (both 
for albedo and for the penetration depth of the laser 
beam). Good quantitative agreement between the re- 
sults of our computer simulation in the active stage 
of the experiment with laboratory data can be treated 
as an additional proof that the required temperature 
distribution inside the sample has been calculated ac- 
curately. 

5. Results 

With the procedure described in the previous Sec- 
tions, we derived the thermal conductivity of the 
three diff'erent dust sample types described in Sect.H 
all consisting of the same 1.5yum-sized Si02 dust. 
For all dust samples we recorded the temperature dis- 
tribution of the sample surfaces by the IR camera for 
ten diff'erent laser intensities. Exemplary, we present 
the results of two compressed dust samples, three 



sieved dust samples, each produced by a mesh with 
a diff'erent aperture size, and two dust samples made 
by the RBD method. In Table[2l each numbered sam- 
ple is listed along with the respective volume filling 
factor (f) and the used laser power P laser- 

Table 2: Parameters of the three different dust sample types: 
compressed (C), sieved (S), and produced by the RBD method 
(RBD). Here, amesh denotes the mesh aperture, the volume 
filling factor of the bulk sample, and Piaser the laser power. 



no. 


sample 





Piaser [W] 


CI 


compressed 


0.54 


1.9 


C2a 


compressed 


0.50 


1.0 


C2b 


compressed 


0.50 


1.9 


SI 


sieved, 


0.29 


1.0 




cimesh = 0.50 mm 






S2 


sieved, 


0.24 


1.0 




cimesh = 0.25 mm 






S3 


sieved. 


0.16 


1.0 




cimesh = 0.15 mm 






RBDl 


RBD 


0.15 


1.9 


RBD2 


RBD 


0.15 


1.9 



As an example of the result of the numerical 
model, Fig. [3] shows the temperature distribution in 
a vertical cross section through sample SI at the end 
of the heating phase. It is readily visible that the heat 
wave penetrates to depths of about 1 mm, which is 
much less than the thickness of the samples. The ra- 
dius of the heated zone exceeds its depth by about a 
factor of three, due to the diameter of the laser beam. 

Fig. |4] depicts the measured temporal and spatial 
surface temperature evolution as a contour plot for 
sample SI (top) in comparison to the respective sim- 
ulated result (bottom). Both diagrams were made by 
averaging the two-dimensional surface temperature 
data in 1 pixel- wide annuli around the center of the 
laser spot. The heating phase during the first 9.5 s 
can clearly be distinguished from the cooling phase 
thereafter. Apart from minor deviations, both fig- 
ures demonstrate the good compliance of measure- 
ment and simulation. The deviations in the active 
phase for distances < 0.5 mm of the numerical con- 
tour plot correspond to the non-gaussian shape of the 
laser profile. 
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Figure 3: Vertical cross section through the temperature distri- 
bution inside the modeled sample S 1 after a heating duration of 
9.5 s. The isotherms are plotted with a division of 10 K. 

In Figs. [51 [6l and [71 the temperature evolution at 
a distance of d = 1 mm from the axis of the laser 
beam is shown for all eight samples listed in Ta- 
ble [2 to visualize and analyze the congruence of the 
experimental (marked as solid lines) and numerical 
(marked as dashed lines) data in detail. These tem- 
perature curves correspond to a vertical cut through 
Fig. [4] at a distance to the center of 1 mm. The 
dashed lines in the figures show the best fitting simu- 
lated temperature curves, achieved by systematically 
varying the value of the thermal conductivity for the 
cooling phase. The other two free parameters, i.e. 
the albedo of the samples and the extinction depth 
of the laser beam, were only slightly varied during 
the heating phase with respect to the values estimated 
by experimental measurements (see Sect. 14.11) . The 
experimentally derived temperature curves show at 
few positions temperature jumps of up to several de- 
grees K. This eff'ect arises from the self-calibration 
of the IR camera to correct a time and temperature 
dependent internal temperature drift. The errors of 
the experimental temperature data, displayed as gray 
shaded areas, are defined by the absolute accuracy of 
the IR camera of ± 2 K and by the size of the tem- 
perature calibration steps, depending on the value 
which is dominating. Lying centered within the er- 
ror intervals of the measured temperature curves, the 
simulated temperature curves indicate a very good 
agreement with the experimental data. For com- 
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Figure 4: Comparison between experimental (top) and numer- 
ical (bottom) results for the temporal and spatial temperature 
evolution of sample S 1 . The contours display the heating (0 . . . 
9.5 s) and the cooling phase (9.5 ... 40 s) of the sample. The 
temperature intervals of the isotherms correspond to 10 K. 

parison, the dotted line for sample C2b in Fig. \5\ 
shows the model results for the measured albedo of 
0.983 and the maximum allowed penetration depth 
of 0.1mm. The resulting thermal conductivity of 
k = 0.078 Wm"^ K"Hs obviously much too high to 
fit both the heating and the cooling curve. 

To show the sensitivity of the simulated temper- 
ature curve on the thermal conductivity as a free 
parameter in the model, the middle Fig. [6] displays 
two additional (dotted) curves with a thermal con- 
ductivity value of 0.0052 Wm"^ K"^ lying above and 
0.0021 Wm"^ K"^ lying below the best compatible 
curve with k = 0.0036 Wm'^ K'^ 

The agreement between of the experimental and 
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Figure 5: Temporal temperature evolution of the compressed 
dust samples CI, C2a, and C2b, showing the heating and cool- 
ing phase for a given distance of 1 mm to the center of the laser 
beam. The solid lines correspond to the results of the experi- 
mental measurements and the dashed lines represent the model 
calculations. The gray shaded area denotes the errors of the ex- 
perimental temperature data. The dotted line for sample C2b 
shows the best-fitting model for the measured albedo of 0.983 
and the maximum allowed penetration depth of 0.1 mm. 



Figure 6: Temporal temperature evolution of the sieved dust 
samples SI, S2, and S3, showing the heating and cooling phase 
for a given distance of 1 mm to the center of the laser beam. 
The solid lines correspond to the results of the experimental 
measurements and the dashed lines represent the model calcu- 
lations. The gray shaded area denotes the errors of the exper- 
imental temperature data. Additionally, the dotted lines in the 
diagram of sample S2 illustrate the numerical results achieved 
by a varied thermal conductivity value of ± 43% of the best fit- 
ting curve. 
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Figure 7: Temporal temperature evolution of the dust samples 
RBDl and RBD2, produced by the RBD method, showing the 
heating and cooling phase for a given distance of 1 mm to the 
center of the laser beam. The solid lines correspond to the re- 
sults of the experimental measurements and the dashed lines 
represent the model calculations. The gray shaded area denotes 
the errors of the experimental temperature data. 



modeled temperature curves in Figs.[5l[6l and [7] pro- 
vides on the one hand a quantitative proof of the va- 
hdity of the model, but serves on the other hand for 
a qualitative analysis of the thermal conductivity of 
each sample. The shape of each individual tempera- 
ture curve comprises some characteristics that argue 
independently for a high or low thermal conductiv- 
ity of the sample. The following quantities suggest a 
low thermal conductivity in comparison to a sample 
with a higher thermal conductivity: 1.) The slope of 
the temperature curve at the start of the heating phase 
is steep, due to the inability of the sample to carry 
away the absorbed energy, 2.) the temperature max- 



imum has a higher value for a fixed energy flux, and 
3.) the temperature curve has a shallower slope at 
the start of the cooling phase for a given temperature 
increase during the heating phase, due to the same 
argument as in 1). If these characteristics are not ful- 
filled, a rather high thermal conductivity value is in- 
dicated. As the temperature maximum is mainly de- 
fined by the absorbed laser power, one should mind 
to compare the maximum values only between sam- 
ples heated with the same laser intensity. 

The comparison between the experimental and nu- 
merical derived data of sample SI in Fig. |4] and the 
temporal temperature evolution of sample S2 (see 
Fig. [8]), depicted for three distances to the center of 
the laser beam (d = 1, 1.5, 2 mm), represent exem- 
plary the stability of the simulated temperature dis- 
tribution in temporal as well as in spatial direction. 
Mind that the thermal conductivity of sample S2 as 
shown in Fig. [8] was derived at a distance of 1 mm 
and still describes the temperature evolution at larger 
distances quite satisfactorily. 
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Figure 8: Temporal temperature evolution of the sieved dust 
sample S2 for three distances (1, 1.5 and 2 mm) to the center 
of the laser beam. The solid lines correspond to the results of 
the experimental measurements and the dashed lines represent 
the model calculations. Mind that the thermal conductivity of 
sample S2 was derived at a distance of 1 mm and still describes 
the temperature evolution at larger distances quite satisfactorily. 

In Table [3l the results for the thermal conductivity 
values of the eight samples are listed, following from 
the best fitting simulated data. As initially expected, 
all samples with relatively low volume filling factors 
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of ^ 0.15 have a thermal conductivity in the or- 
der of 10"^ W m"^ K"\ whereas the samples with the 
highest (p values ^ 0.50) have an order of magni- 
tude higher values for the thermal conductivity. 

Table 3: Values of the model calculations for the thermal con- 
ductivity the albedo, and the extinction length A, which pro- 
vide the best agreement with the experimentally derived tem- 
perature evolution of the eight dust samples in comparison to 
their volume filling factors (p. 



sample 




k [W K-i] 


albedo 


A [yum] 


CI 


0.54 


0.021 


0.9960 


50 


C2a 


0.50 


0.021 


0.9935 


50 


C2b 


0.50 


0.016 


0.9945 


50 


SI 


0.29 


0.0078 


0.965 


50 


S2 


0.24 


0.0036 


0.980 


100 


S3 


0.16 


0.0016 


0.972 


100 


RBDl 


0.15 


0.0026 


0.992 


150 


RBD2 


0.15 


0.0026 


0.990 


100 



An estimate of the accuracy with which the ther- 
mal conductivity is measured in our work is a non- 
trivial task. We can confidently discuss the accuracy 
of the temperature fields measured in the laboratory 
experiments and the accuracy of the computer simu- 
lations. However, both do not completely determine 
how accurately we estimated the thermal conductiv- 
ity. The development of a statistical algorithm of ac- 
curacy estimation for the proposed technique is be- 
yond the scope of this work. At present, we use the 
empirical evaluation approach described below. 

Let us first quantify the allowed range of the model 
parameters and their relative errors. We begin with 
the estimates of the albedo. As previously men- 
tioned, the albedo of the compressed sample was 
experimentally determined. The measured value is 
only 2% smaller than our model value (see Table [3]). 
However, it is important to note that it is the albedo 
(i.e., the intensity of the scattered radiation) that was 
measured, whereas the absorbed energy, directly pro- 
portional to (1 -albedo) plays a decisive role in our 
computer modeling of the heat transfer. It is obvious 
that for such a high albedo, even minor changes dra- 
matically alter the amount of absorbed energy and, 
consequently, the temperature evolution of the sam- 



ple (under the condition that the other model param- 
eters are fixed). The second free parameter of the 
model is the attenuation factor of the laser radiation 
(e-folding scale). Accurate measurements of the at- 
tenuation of radiation in porous media is a very diffi- 
cult task. However, we can get some quantitative es- 
timates based on the experimental evaluation of the 
transparency of our samples. Our estimates show 
that for the compressed samples, even for a sample 
thickness of a small fraction of a mm, no more than 
a few percent of the absorbed energy passes through. 
This gives us an upper limit of the attenuation factor 
of ~ 0.1 mm. 

Firstly, we tested the possibility to find a better 
agreement between experiment and computer model 
using the measured values of the albedo. As the 
amount of the absorbed energy is much larger in this 
case, one has to increase the value of the eff'ective 
thermal conductivity in the model in order to fit the 
experimental data. Based on the law of conserva- 
tion of energy, one could expect that the new value 
of the coefficient of thermal conductivity should be 
roughly three times larger than the value presented 
in the Table [3l Our attempts to find a good match be- 
tween measured and calculated temperature curves 
show that for all physically reasonable model param- 
eters (attenuation factor < 0.1 mm and reduction fac- 
tor h < 0.05) the standard deviation is significantly 
greater if we use the measured albedo rather than the 
value obtained for the "optimal" set of parameters. 

Thus, we are confident that the values of the ther- 
mal conductivity shown in Table [3] are in reason- 
able agreement with the "true" thermal conductiv- 
ity of the samples. As stated above, absolute er- 
rors are hard to give; however, from the variation of 
the model thermal conductivities for identical sam- 
ples, we guess that the relative error of our derived 
thermal-conductivity values should be < 30%. 

If we plot the thermal conductivity values 
from Table [3] and the value for fused silica of 
1.4 Wm"^ K"^as a function of the bulk volume fill- 
ing factor, we get the behavior depicted in Fig. [9l 
Together with the value of silica glass, the thermal 
conductivity seems to be exponentially dependent on 
the volume filling factor of the sample with the cor- 
responding relation = 0.000514^^-^^^ Wm"^ K'^ 
Further measurements with diff'erent materials and 
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Figure 9: Thermal conductivity of fused silica and all dust 
samples studied in this work plotted versus their volume fill- 
ing factor. The dashed line displays the fitted exponential law 
k = 0.000514^^-91'^ Wm-i ^-i^ 



porosities should close the gap for higher volume fill- 
ing factors and verify this exponential correlation. 

6. Discussion 

We developed a measuring technique to non- 
invasively investigate the thermal conductivity of 
very porous and fragile dust samples. A combination 
of experimental measurements and numerical sim- 
ulations, precisely emulating the laboratory experi- 
ments, provides the determination of the heat con- 
ductivity of the sample. 

In this work, we determined the heat conductiv- 
ity of samples consisting of the same bulk material 
(monodisperse Si02 spheres with 1.5yum diameter) 
but with diff'erent inner structures, resulting from dif- 
ferent manufacturing techniques (RBD, sifting, com- 
pression). Our measurements and the according sim- 
ulations could resolve the thermal conductivity for 
sample volume filling factors of (p = 0.15 to 
^ = 0.54, showing by trend an exponential correla- 
tion between the thermal conductivity of the sample 
and its porosity, with conductivity values as low as 
k ~ 10"^ Wm"^ K"^ for the samples with the lowest 
packing density. 

Regarding the heat conductivity as a means of en- 
ergy transport, the heat conduction is determined by 
the contact area between adjacent single particles of 
the bulk material. As the dust used in this study has a 



very narrow size distribution, we assume equal-sized 
contact areas for all contact points. Thus, the num- 
ber of contacting grains (coordination number) influ- 
ences the efficiency of heat transport within the mate- 
rial. The coordination number of a material depends 
on the microscopic and macroscopic homogeneity of 
the structure. In the case of our samples, we assume 
a homogenous distribution of single particles within 
the compressed samples and the samples made by the 
RBD method. In these cases, the assumption that the 
thermal conductivity depends directly on the coordi- 
nation number and, thus, the porosity seems to be 
valid. For samples with a more complex inner struc- 
ture, like the sieved dust samples, the correlation 
of thermal conductivity, porosity and coordination 
number is more complex. The sieved dust samples 
consist of a combination of porosities: the volume 
filling factor of the sifted dust aggregates (^agg) and 
the packing density of the sifted aggregates inside the 
repository (0««). The heat conductivity within a sin- 
gle sieved dust aggregate diff'ers from the heat con- 
duction between these aggregates, due to more con- 
tact points in between the monomers inside individ- 
ual sub-mm-sized aggregates. Fig. [TOl illustrates the 
distinction of the heat conductivity dependent on the 
inner structure of the dust sample. The sketch dis- 
plays stepwise the heat conduction between adjacent 
monomers. Each time step comprises the heat trans- 
port through one neighboring monomer layer. By 
comparison, the compressed dust sample shows the 
most particles, which were involved in the heat trans- 
port, as they obviously have the most contact points 
between monomers. If we consider the sieved and 
the RBD dust sample, both show approximately the 
same amount of heated single particles, but the RBD 
sample slightly more. If one supposes that the ther- 
mal conductivity only depends on the porosity, one 
would expect the contrary: that the sieved dust sam- 
ple with a lower porosity compared to the RBD sam- 
ple would transfer the heat better. The sketch and the 
measured results of sample S3 compared to the sam- 
ples RBDl and RBD2 indicate that not the porosity 
but the number of contacting particles has the ma- 
jor influence on the thermal conductivity. From time 
step to time step the RBD dust sample has uniformly 
only a few adjacent particles which transfer the heat. 
In contrast, the sieved dust sample has local regions 
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Figure 10: Sketch to illustrate the stepwise heat conduction in- 
side the three differently structured sample types. Left column: 
compressed dust sample; middle column: sieved dust sample; 
right column: dust sample produced by the RBD method. 



(within the sifted dust aggregates) where many con- 
tacting particles transfer the heat in one time step 
whereas between these aggregates only rare contact 
points exist, resulting in a less efficient heat conduc- 
tion compared to the RBD sample. 

We performed the experiments described in Sect. [5] 
under high-vacuum conditions so that the thermal 
conductivity of the residual gas was negligible. It 
is not a priori clear whether this is also true for solar- 
nebula conditions. The thermal conductivity of a 
diatomic gas with pressure p and temperature T is 
given by 



5 pul 
^ 6 



(5) 



with u and / being the average molecular veloc- 
ity and the inter-particle distance within the dust 



sample, respectively. For our dust samples we ex- 
pect / = 1...10yum. To estimate the contribu- 
tion of the gas to the thern ial conductivity we use 
the solar nebu la mode l s bv I Weidenschilling ( 1977 ) 
fMMS N) and IPeschl (|2007Tk Following the val- 



ues of IZsom et al.l (120101) for the midplane den- 
sity of a protoplanetary disk at 1 AU of Pmmsn = 
1.4 X 10"^gcm"^ and poesch = 2.7 x 10"^gcm"^ 
for T = 200 K, we get for the thermal conductiv- 
ity of the gas kg^MssN = 7.0 x 10"^ Wm"^ K"^ and 
kg,Desch = 1.3 X 10""^ Wm"^K"^ for / = I fim and 
kg,MssN = 7.0 X 10"^ Wm"^ K"^ and kg^oesch = 1-3 x 
10"^Wm"^K"^ for / = lOyum, respectively. To 
evaluate the influence of the gas conductivity on the 
total thermal c onductivit y, we applied the model b\ 



total tnermal c onductivit y, we appliea tne model by 
Russelll (1 1935b for which Icheng and VachonI (Il970[) 



found a good agreement with measurements. For our 
compact samples, the increase of the total conductiv- 
ity due to the gas filled pores is less than 4% for the 
whole range of solar nebula gas densities. The in- 
fluence of the gas on the total conductivity becomes 
slightly more important for those samples with high 
porosity and low thermal conductivity. The thermal 
conductivity of sample S3 is increased by 75% if 
we assume a thermal conductivity of the solar neb- 
ula gas of kg^Desch = 1-3 X 10"^ Wm"^ K"^; all other 
gas densities result in a thermal conductivity increase 
of less than 10%. Thus, we can conclude that the 
influence of the gas on the thermal conductivity is 
negligibly small for all cases but the highest gas den- 
sities and the lowest volume filling factors, and even 
then only a small correction is required. The role of 
gas may be m ore important in the case of porous ice. 
For instance, iBar-Nun and Lauferl (|2003b measured 
for amorphous "cometary" water ice samples an in- 
crease of the thermal conductivity by a factor of 3-4 
due to the eff'ect of internal trapped gas. 

The measurement of the thermal conductivity and 
its dependence on the porosity has significant impli- 
cations for the thermal evolution of meteorite parent 
bodies in the early Solar System. Presently, it seems 
clear that meteorite parent bodies accreted within 
less than a few Myrs after the first solid objects - 
Ca,A l rich i nclusions - formed (e.g. Trieloff (|2009b : 
Scottl(l2006b ). Hence, short-lived nuclei, mainly Al 
and ^^Fe, were eff'ective energy sources to heat early 
formed planetesimals. Classical thermal models 
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(IMiyamoto et allll98ll : lGhosh and McSweenl . ll998b 
apply thermal conductivities measured on mete- 
orites, which are already solidified rocks with only 
small degrees of porosity. However, meteorite parent 
bodies started with likely higher porosity, so heating 
eff'ects can be expected for comparatively small plan- 
etesimals. This potentially implies significant heat- 
ing eff'ects of smaller bodies during their growth to 
asteroid sizes, which would imply drastically diff'er- 
ent thermal histories and has to be taken into account 
when modeling or describing the history of planetary 
building blocks in the early Solar System. Hence, 
our measured thermal conductivities should be im- 
plemented into models that describe the early heat- 
ing of meteorite parent bodies, i.e. asteroids. For 
comets, which retained highly volatile species like 
CO, CO2, and H20-ice particles, as recently discov- 
ered in Comet Hartley-2 (A'Hearn et al., in prepa- 
ration), heating eff'ects were obviously very minor, 
so that they retained their originally porous structure 
that likely was also a property of meteorite parent 
bodies before heating by short-lived nuclide decay 
energy. The pristinity of comets in terms of porous 
structure and low thermal conductivity is demon- 
strated by the agree ment of the thermal conductiv- 
ity values derived by lOroussin et al.l OOOVb for comet 
9P/Tempel 1 {k < 1.8x10"^ . . . 3.6x10"^ Wm'^ R-^) 
and the thermal conductivity of our low-density sam- 
ples (fc= 1.6x 10-^..2.1 X 10-2 Wm-i R-^). 
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